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^ ! 1 Introduction 

Spatial patterns of gene expression are central to the development of multi-cellular organisms. Most math- 
• ematical studies of pattern formation investigate diffusion-driven instability, which is a mechanism that 

Q ! amplifies spatial inhomogeneities in a class of reaction-diffusion systems (see, e.g., HI). However, many 

ch ! patterning events in multi-cellular organisms rely on cell-to-cell contact signaling, such as the Notch path- 

way 0, and do not involve diffusible proteins for intercellular communication. A particularly interesting 
phenomenon in this form of communication is lateral inhibition whereby a cell that adopts a particular fate 
inhibits its immediate neighbors from doing likewise Q, thus leading to 'fine-grained' patterns. There is 
increasing interest in understanding the Notch signaling circuitry in mammalian cells that leads to such lat- 
eral inhibition 131|5l. Recent studies showed that a lateral inhibition pathway also functions in E. Coli, and 
t^- , enables the bacteria to inhibit the growth of other E. Coli strains in direct contact Q. 
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Dynamical models are of great interest for understanding the circuit topologies involved in lateral inhi- 
bition and for predicting the associated patterns. Several simplified models have been employed for Notch 
signalling pathways in 01 and Q. The objective of this paper is to present an abstract dynamical model that 
captures the essential features of lateral inhibition and to demonstrate with dynamical systems techniques 
that these features indeed lead to patterning. Although this model is not meant specifically for Notch signal- 
ing, it encompasses as special cases the lateral inhibition model in Q as well as a slightly modified version 
of the one in Q. 

Our model treats the evolution of concentrations in each cell as an input-output system, where the inputs 
represent the influence of adjacent cells and the outputs correspond to the concentrations of the species that 
interact with adjacent cells. The input-output models for the cells are then interconnected according to an 
undirected graph where the nodes represent the cells, and the presence of a link between two nodes means 
that the corresponding cells are in contact. The main assumption on the input-output model is that each 
constant input yields a unique and globally asymptotically stable steady-state, and that the value of the output 
at this steady-state is a decreasing function of the input. This decreasing property captures the inhibition 
of the cell function by its neighbors. The model allows multiple inputs and outputs, and is restricted by a 
monotonicity assumption, following the definition of monotonicity for dynamical systems with inputs and 
outputs Q. 
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Using this model, we first give an instability condition for the homogeneous steady-state, applicable to 
arbitrary contact graphs. We then focus our attention on bipartite graphs, and demonstrate the emergence 
of a "checkerboard" pattern, exhibiting alternating high and low values of concentrations in adjacent cells. 
Next, we establish a strong monotonicity property of the interconnected model for bipartite graphs, which 
implies that almost every bounded solution (except for a measure-zero set of initial conditions) converges 
to a steady-state (SJI3- A graph is bipartite if and only if it contains no odd-length cycles, and Cartesian 
products of bipartite graphs are also bipartite ifTOll . Thus, the results of this section are applicable, among 
others, to grid graphs (one dimensional path graphs and their Cartesian products in higher dimensions) 
which are appropriate for representing arrays of cells. 



2 Lateral Inhibition Model and Preliminaries 

We let Q be an undirected, connected graph where the nodes represent the cells, and the presence of a link 
between two nodes means that the corresponding cells are in contact. In preparation for the dynamical model 
studied below, we let N denote the number of cells and define the matrix P e R NxN : 

( di~ { if nodes i and j are adjacent, 
1 otherwise, 

where dt denotes the degree of node i. It follows that P is a nonnegative row-stochastic matrix, that is: 

PI = 1 (2) 

where 1 denotes the vector of ones. The matrix P is identical to the probability transition matrix for a 
random walk on the graph Q. The properties summarized below therefore follow from standard results for 
random walks (see, e.g., ifTTI ): 

Lemma 1. P possesses real eigenvalues A^ < ■■■ < A\ all of which lie in the interval [—1,1], and corre- 
sponding real, linearly independent eigenvectors v,-, i = l,-- ,N. In particular, A\ = 1, and v\ - 1 is a 
corresponding eigenvector. IfQ is bipartite, then A^——\, and an eigenvector v^ is such that the entries are 
either 1 or — 1, and two entries corresponding to adjacent nodes have opposite signs. 

Let i - 1, • • • , N denote the cells, and consider the dynamical model: 

x l = f(x i ,u i ) y = h(x') (3) 

where x' e 2£ c R" is a vector describing the state of reagent concentrations in cell i, u' e ^ c R m describes 
the 'input' from adjacent cells, and y' e <3/ c W" describes the 'output' that serves as an input to adjacent 
cells. In particular, 

U = (P®I m )Y (4) 

where P is as defined in £[]), U := [u lT ■ --u nT ] t and Y := [y lT ■ ■ -y NT ] T . If follows from CD) that the input u' 
is the average of the outputs y k over all neighbors k of cell i. Thus, we henceforth take the input and output 
spaces to be identical: % = & . 

We assume that /(■, ■) and h(-) are continuously differentiable and further satisfy the following property: 

Assumption 1. For each constant input u*, system (0) has a globally asymptotically stable steady-state 
x* := S(u*) with the additional property that: 



j df(x,u) 



(x,u)=(x* ,u*) 



tO. (5) 
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The map S : % — > X and, therefore, the map T : — > defined by: 

T(-):=h(S(-)), (6) 

are continuously dijferentiable. 

Following the terminology in Q, we will refer to S(-) as the input-state characteristic, and to T(-) as 
the input-output characteristic. Our next assumption is that (f3]) is a monotone system in the sense of 0, as 
defined below. According to the classical definition for systems without inputs and outputs J9l, a monotone 
system is one that preserves a partial ordering of the initial conditions as the solutions evolve. The partial 
ordering is denned with respect to a positivity cone K in the Eucledean space that is closed, convex, pointed 
(K n (-K) = {0}), and has nonempty interior. Given such a cone, x < x means x-xeK,x<x means x < x 
and x + x, and x «: x. means that x. - x is in the interior of K. The system x = f(x) is then defined to be 
monotone if two solutions x(t) and x(t) starting with the order x(0) < x(0) maintain x(t) < x(t) for alfl t > 0. 
The more restrictive notion of strong monotonicity stipulates that x(0) -< x(0) implies x{t) <sc x{t) for all t > 0. 
The monotonicity concept was extended to systems with inputs and outputs in Q : 

Definition 1. Given positivity cones K U ,K Y ,K X for the input, output, and state spaces, the system x = 
f(x, u), y - h(x) is said to be monotone if x(0) < x(0) and u(t) < u(t) for all t > imply that the resulting 
solutions satisfy x{t) < x(t)for all t>0, and the output map is such that x < x implies h(x) < h(x). 

Assumption 2. The system (0) is monotone with respect to K u = R™ Q , K Y = -K u , and K x = K, where K is 
some positivity cone in K". 

As observed in Remark V.2], monotonicity implies that the input-state and input-output characteris- 
tics are nondecreasing with respect to the same ordering; that is, u < u with respect to K u implies S(u) <S (u) 
with respect to K x and T(u) < T(u) with respect to K Y . Since K Y = -K u in Assumption |2j we conclude that 
T(-) is nonincreasing with respect to the standard order induced by K u = R'" Q . This nonincreasing property 
means that, if two cells are in contact, an increase in the output value of one has the opposite effect on the 
other, which is why ©-dU) is referred to as a "lateral inhibition" model. We note from the nonincreasing 
property of T(-) that: 

T\u) := — K — (7) 
ou 

is a nonpositive matrix in R mXm ) and denote its spectral radius as: 

p{T\u)). (8) 

We conclude this section by quoting lemmas that will be used in the sequel. Lemmas |2] and [3] are from lfT2l : 

Lemma 2. Given the system x - f(x,u), y = h(x) with continuously differentiable /(•,•) and h{-), the lin- 
earization x -Ax + Bu, y = Cx about a point (x*,u*) satisfying f(x*,u*) = is also monotone with respect 
to the same positivity cones. 

Lemma 3. The linear system x - Ax + Bu, y = Cx is monotone if and only if: 

1) xe K x implies Ax e K x , 

2) ueK u implies Bu e K x , 

3) xe K x implies Cx e K Y . 



'Here, "for all f" is understood as "for all times / that belong to the common domain of existence of the two solutions." 
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The following lemma, proven in lfT2l for single-input, single-output systems and extended in [13] to the 
multivariable case, determines stability of a positive feedback system based on the 'dc gain' of the open-loop 
system: 

Lemma 4. Suppose the linear system x - Ax + Bu, y = Cx is monotone with respect to cones K U ,K Y ,K X 
such that K u = K Y and A is Hurwitz. If —(I + CA~ l B) is Hurwitz, then so is A + BC. If—(I + CA~ l B) has an 
eigenvalue with a positive real part, then so does A + BC. 

In the special case of single-input, single-output systems, the stability condition above amounts to checking 
whether the dc gain -CA~ l B is greater or smaller than one. In the multi-input, multi-output case, this 
condition is equivalent to inspecting whether the spectral radius of the dc gain matrix is greater or smaller 
than one. 

The following test from QQl is useful for certifying monotonicity with respect to orthant cones: 

Lemma 5. Consider the system x = f(x, u), y - h(x), x e f c P, u e % c W, y e *3f cP, where 
the interiors of X and % are convex, and /(•>•) and h(-) are continuously differentiable. If there exist 
€i,--- ,€„,Su--- ,5 m ,pu--- ,p p e {0,1} such that: 

Of. 

ox k 

(_l)^+<5*|/j( X)M )>0 \/ xe ^,\/ue^,\/j,k (10) 
du k 

(-l) e ' + " l ^(x,B)>0 Vxe3T,\fj,k, (11) 

dXj 

then the system is monotone with respect to the positivity cones K u = {u 6 R" ! | (— V) J Uj > 0], K x = {x e 
R" | {-lf ] xj > 0}, K Y -{yeW \ (-Vf'yj > 0}. 

3 Instability of the Homogeneous Steady-State 

Note that system ©-(HI) admits spatially homogeneous solutions of the form x'(t) = x(t), i= 1, • • • ,N, where 
x(t) satisfies: 

x = f(x,h(x)). (12) 

In particular, if the map T(-) has a fixed point: 

u* = 7\u*), (13) 

then (fT2l admits the steady-state: 

x*=5(u*)- (14) 

For single-input, single-output systems with % = R>o, the nonincreasing property of the map T : R>o — > R>o 
indeed guarantees a unique fixed point u* in (fT3l . 

The "lumped model" (fT2l describes the dynamics of the A^n-dimensional system © reduced to the n- 
dimensional invariant subspace where the solutions are spatially homogeneous. Thus, the steady-state x* of 
the lumped model defines the homogeneous steady-state x l - X*, i - 1, • • • ,N, for the full system ©-(HJ). As a 
starting point for the analysis of pattern formation, we now give an instability criterion for the homogeneous 
steady-state: 
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Theorem 1. Consider the system (0)-@ and suppose Assumptions\I\and\2\hold. Let An denote the smallest 
eigenvalue ofP as in Lemma\J} and let u*, x* be as in HI J I ), < \14i . If: 



A NP (T'(u*)) < -1, (15) 
then the homogeneous steady-state x' -x*, i - I, - ■ ■ ,N, is unstable. 

Proof: Let X := [x iT ■ • ■ x nT ] t , and note that the linearization of ©-(|4]) about the homogeneous steady-state 
[x* T 



,x* T ] T gives the Jacobian matrix: 



where: 



A := 



df(x, u) 



dx 

We recall from Lemma[T]that 



I N ®A + P®{BC) 
df(x, u) 



B 



(x,m)=(x*,u*) 



y- l PV = A:-- 



du 



C:-- 



dh(x) 



(x,k)=(x*,u*) 



dx 



An 



(16) 
(17) 

(18) 

(19) 

(20) 
(21) 

then (120\) has a positive eigenvalue. 

The theorem follows from this claim because, if (fT5T ) holds, then (l20l has a positive eigenvalue for 
k - N, which implies instability. To prove the claim, we note from Assumption |2] and Lemma |2] that the 
linear system x = Ax + Bu, y = Cx is monotone with respect to K u = R'" , K Y = -K , and K x = K. We write 
A + AkBC -A + BCk where Cu := A^C and note that (|2TI ) implies A& < 0. Thus, the system x = Ax + Bu, y = C^x 
is monotone with K u = K Y . In addition, Assumptions Q] and |2] imply that A is Hurwitz, as can be deduced 
from Ifl2l Lemma 6.5]. Thus, it follows from the second statement of Lemma @] that if -(/ + C^A -1 B) has a 
positive eigenvalue, then so does d20l) . The remaining task is thus to prove that 



where V = [vi ■ ■ • vjy], and apply the following similarity transformation to (1161 ): 

(V- 1 ® I n )[I N ® A + P ® (BC)](V®I n ) = I N ®A + A® (BC). 
This matrix is block-diagonal, with the kth diagonal block given by: 

A + A k BC. 

Claim: If 

A kP {T\u*))<-\, 



-(1 + C k A' l B) = -I- A k CA~ l B 
has a positive eigenvalue. To this end, we first show that 



Since 



differentiation gives: 



df(x, u) 



dx 



T'(u*) - -CA~ l B. 



f(S(u),u) = 0, 



dS (u) df(x, u) 
x=s(u) du du 



0. 



x=S («) 



(22) 

(23) 
(24) 
(25) 
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Next, it follows from the definition © that 

dS (u) 



T (u) 



(26) 

x=S(u) OU 



dx 

Combining (T25T ) and (l26l >. and substituting (fTTT ). we verify (T23T ). Substituting (1231) . we then rewrite (1221 as 

-I + A k T'(u*), (27) 

and conclude that it indeed has a positive eigenvalue, because A k < implies that A^T'iu*) is a nonnegative 
matrix and (|2TI ) implies that its spectral radius exceeds one. Since the spectral radius is an eigenvalue for 
nonnegative matrices (see, e.g., lfT5l ). the conclusion follows. □ 

The eigenvectors v k of P used in the similarity transformation ( fT9l ) may be interpreted as the spatial 
modes of the system. Thus, the stability properties of the matrix (1201 ) for each k determines whether the 
corresponding mode decays or grows in time. Since the spectral radius is nonnegative and k — 1,- • • ,N, 
are in decreasing order, whenever the instability criterion (|2TT) holds for a particular mode k, it also holds for 
higher values of k. Because larger wavenumbers k imply higher spatial frequency content in vu, we conclude 
that the instability condition above sets the stage for the formation of high-frequency spatial patterns. 

4 Patterning in Bipartite Graphs 
4.1 Emergence of Checkerboard Patterns 

For bipartite graphs, where A^ = -1 as stated in Lemma [T] the instability condition in Theorem Q] is: 

p(T'(u*)) > 1. (28) 

This condition indicates the growth of the highest spatial-frequency mode which exhibits opposite signs 
for adjacent nodes. Thus, concentrations in adjacent nodes move in opposite directions in the vicinity of the 
homogeneous steady-state. We now show that, if the map 

r 2 (-) := 7X7X-)) (29) 

has two fixed points ui + m other than u*, satisfying: 

ui = 7Xu 2 ), u 2 = 7Xui), (30) 

then the system ©-(lU) has an inhomogeneous steady-state with two sets of concentrations, each assigned 
to one of two adjacent cells. We will refer to this steady-state as a "checkerboard" pattern, since adjacent 
cells adopt distinct states. Although this term may be associated with cells arranged as a grid graph in two 
dimensional space, we will use it broadly for any spatial arrangement that forms a bipartite graph. 

Proposition 1. Let Qbe a bipartite graph and let the sets I c { 1 , • • • ,N] and I' = { 1 , • • • ,N} — I be such that 
no two nodes in the same set are adjacent. If there exist ui e ^ and 112 £ % , Ui ^ 112, satisfying A30]) . then 

x 1 = S (ui), i £ I, x' =S(u 2 ), i e (31) 

and 

x' - 5(u 2 ), i e I, x 1 = S(m), i e I', (32) 
are steady -states for system ©-(E]). 
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Proof: To show that (I3TT ) is a steady-state, we note that, if i e I, then / = T(u\) and, if i e I', then y' = 7\ii2). 
From ©, the input u l to a node in I is T(u2) because all neighbors of this node belong to I' . Likewise, 
the input u' to a node in I' is T(u\) because all neighbors of this node belong to I. Since T(u2) = Ui and 
7\ui) = U2, we conclude that (I3TI ) is indeed a steady-state, and identical arguments apply to (l32l ). □ 



Theorem 2. Consider the system (0)-dU) and suppose Assumptions\J\and^\ and the hypotheses of Proposi- 
tions/told. If, in addition, 

P (r(ui)r'(u 2 )) < l, 03) 

then the steady-states (13 li and (132\) are asymptotically stable. 

Before giving the proof, we note that (1301 ) corresponds to a period-two orbit of the discrete-time system: 

u(t+l) = T(u(t)), (34) 

and (1331 implies the asymptotic stability of this orbit. Likewise, (1281) indicates instability of the fixed point 
u* for this discrete-time system. Thus, an interesting duality exists between (l34l and the spatially-distributed 
system ©-(HJ) defined on a bipartite graph: A bifurcation from a stable fixed point to a stable period-two 
orbit in (l34l corresponds to the emergence of stable checkerboard patterns from a homogeneous steady-state 
in ©-©. 

In the single-input, single-output case with & = R>o, where T : R>o — > R>o is a nonincreasing function 
by Assumption^ condition (1281) indeed implies the existence of a period-two orbit d30l) . To see this, assume 
to the contrary that u* is the unique fixed point of T 2 {-). Since 7\-) is continuous and nonincreasing, this 
uniqueness property would imply that u* is a global attractor for all solutions of the difference equation (1341) 
starting in R>o lfl6l Lemma 1.6.5]. This, however, contradicts (1281) . which implies instability of u* for this 
scalar difference equation. 

The argument above does not suggest the uniqueness of the pair (111,112), and multiple pairs satisfying 
(l30l may exist. However, we claim that at least one pair satisfies: 



dT 2 (u) 



du 



dT 2 (u) 



U=M\ 



r'(ui)r'(u 2 ) < 1, (35) 



»=U2 

which is the scalar equivalent of d33l , since T'(ui)T'(u2) is nonnegative. To see this, note from d28l ) that: 

dT 2 (u) 



du 



= r'(u*)7"(u*) > 1 (36) 



and suppose, in contrast to d35l ), that the derivative of T 2 (-) is greater than or equal to one at each of its fixed 
points. This implies that T 2 (u) > u for all u > u*, because T 2 (u) - u has nonnegative slope at zero-crossings 
and, thus, remains nonnegative for u > u*. The inequality T 2 {u) > u implies unbounded growth of T (■) 
which is a contradiction because T(-) is continuous and nonincreasing, thus, bounded. 

Proof of Theorem |H Let Nj and Nj> :=N — Nj denote the cardinalities of the sets I and I', and index the 
cells such that i - 1, • • • , Nj belong to I, and i = Nj + 1, •• ■ ,iV belong to I'. Then the matrix P has the form: 



P = 



Pu 
P21 



(37) 
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where P l2 £ R NlXNl ' , P 2 i e R Nj,xNl . Let X := [x lT ■ ■■x nT ] t , and note that the linearization of ©-© about 
(f3TT > gives the Jacobian matrix: 



where 
A 



df(x, u) 



j •" 



dx 



B 



I Nl ®Ay Pn®(B l C 2 ) 

P 2 i®(S 2 Ci) I Nj ,®A 2 
df(x,u) 



(X,U) = (S (Uj ),U;) 



./ •" 



du 



C, 



dh(x) 



dx 



x=S(Uj) 



1,2. 



(38) 



(39) 



(x,u)=(S (Uj),Uj) 

From the definition £T|), the matrix DP, where D is a diagonal matrix of the node degrees, is symmetric. 
Since D~ l/2 (DP)D~ l/2 = D l/2 PD~ l/2 is also symmetric, we write: 



D l/2p D -l/2 



R 
R T 



(40) 



where R e R N i xN i'' is appropriately defined. Then, we apply the following similarity transformation to (l38l) : 

I Nl ®Ay 



(D 1/2 ®/„) 



P2l®(B 2 Ci) I Nt , Q9A 2 



(D 



-1/2, 



'/„) - 



R T ®(B 2 C l ) h 



^12®(5iC 2 ) 

The structure of d40l is such that it can diagonalized with an orthonormal matrix of the form: 

Q 



Qi Qi Q3 
Qi -Qi Q A 



(41) 



(42) 



which results in: 



R 
R T 



Q = Q 



-A+ 







(43) 



where A + is a diagonal matrix of the strictly positive eigenvalues of P, the columns of <2 3 and Q 4 span 
the null spaces of R T and R, respectively, and the dimensions of the zero diagonal blocks in (l43l) are con- 
sistent with the dimensions of these null spaces (which we denote as nj, and n 4 , respectively). From the 
orthonormality of Q, we get the identities: 



Q\Q\ = Q T 2 Qi = \h 

e[03 = o q t 2 qa = o, 

where r is the dimension of A + . Likewise, equation d43l ) implies: 



7M 



RQi = QiA, 



R 1 Qi = QiA + 



RQ\ = R T Q 3 = 0. 

We now return to the Jacobian matrix (|4TT) and further apply the following similarity transformation: 



(44) 

(45) 
(46) 



(47) 
(48) 



2Q[®I n 






2Q T ? ®I n 



Q\®In 

Ql®I„ 



I Nl ®A { R®{B X C 2 ) 
R T ®(B 2 Ci) I Nr ,®Ai 



Q\®h Q 3 ® I n 

Q 2 ®I n Q A ®I n 



(49) 
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where the leftmost matrix is the inverse of the rightmost matrix from (I44l- (l46l) . Likewise, using (l44l-(|48l). 
it is not difficult to show that the product d49l equals: 



I r ®Ai 
A + ®(B 2 Ci) 



A + ®(5iC 2 ) 
7 r ®A 2 



'"3 



(50) 



Since Assumptions[T]and|2]imply that A\ and A 2 are Hurwitz Ifl2l Lemma 6.5], stability of (TSUI) is determined 
by the upper left blocks which, upon a similarity transformation with an appropriate permutation matrix, are 
block-diagonalized into r blocks of the form: 



Ai AiB x C 2 
AiB 2 C x A 2 



(51) 



i-l,--- ,r. 

We will now show that (f5Tb is Hurwitz for any Ai e [-1, 1]. Since all eigenvalues of P lie in this interval 
by LemmaCD this will conclude the proof. We do not provide a separate proof for the asymptotic stability of 
(1321 ). as identical arguments apply when the indices 1 and 2 are swapped in (T5Tb . If Ai = 0, (1511 ) is Hurwitz 
because A i and A 2 are Hurwitz. If Ai ± 0, then we apply the similarity transformation: 



/ 


Ai 


AjBiC 2 


I 




Ai 


At 1 ! 


AiB 2 C\ 


A 2 


AJ 




B 2 C X 



A]B X C 2 
A 2 



and rewrite the result as: 



where 



Ax A]B X C 2 " 


, % = 





A 2 




. B 2 



C = [d 0]. 



(52) 



(53) 



(54) 



">0' 



We claim that the linear system defined by the triplet (C, J{, S) is monotone with respect to K u = K 1 
and K x = -KxK where K is as in Assumption^ To see this, first note from Lemma|2]that {C\,A\,B\) and 
(C 2 ,A 2 ,B 2 ) are monotone with respect to the cones specified in Assumption [2 By Lemma [3j this means 
that: 

x e K => Ajx 6 K, ueR™ ^>BjueK, xeK CjxeR™ , 7 = 1,2. (55) 
We now show that the conditions of Lemma[3]hold for (C,Jl,S) with K u = K Y = '. 



1) Suppose x = [x^x^] 7 e —K x K, that is x\ e -K, x 2 e K. Then, 



->() 



,K x = -KxK: 



ftx = 



A ixi +AjB\C 2 x 2 
A 2 x 2 



e -KxK 



(56) 



because, from d55l ), A\X\ e -K, A 2 x 2 e K, C 2 x 2 e and, hence, B\C 2 x 2 e -K. 

2) We want to show that u e W" implies Su e -K x K. From the definition of S in d54l ), Su e -K x K means 
B 2 u £ K. It follows from the second implication in (1531) that u e R" ! indeed implies B 2 u € K. 

3) To prove monotonicity with K Y = R" ! , we need to show that x\ e -K and x 2 e K imply Clx^x 7 ] 7 e R.'" () . 
This is indeed true, since Clx^x 7 ] 7 = C\X\ and, from (1551) . ;q 6 -K implies C\x\ e R^ Q . 
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Having verified the conditions of Lemma |3l we conclude that (C, Si, S) is monotone with respect to 
K u = K Y = K.'" . In addition, the matrix 3\ in d54l ) is Hurwitz, as A\ and A 2 are Hurwitz. Thus, it follows 
from the first statement in Lemma|4]that, if -(I + CJl~ i S) is Hurwitz, then so is (|53T ). Note that 

C$T 1 B=[Ci 0] 

and, from a derivation similar to the one for 423), T'(Uj) = -CjAfB h j = 1,2. Thus, © gives: 

-(/ + C^ _1 S) = -/ + ^r'(ui)r'(u 2 ), (58) 

and (1331 and Aj e [-1, 1] imply that -(l + C^\~ l< B) is indeed Hurwitz. From Lemma|4j this means that (l53l 
and, thus, d5TT > is Hurwitz i = 1, • • • , r, concluding the proof. □ 

4.2 Generic Convergence to Steady-States 

Thus far we have studied local asymptotic stability properties of the steady-states. Strongly monotone sys- 
tems (as defined in the paragraph above Definition []} have been shown to possess a "generic convergence" 
property (8j |9l which means that almost every bounded solution (except for a measure-zero set of initial 
conditions) converges to the set of steady-states. Below we first prove monotonicity of ©-(HJ) in Theorem 
|3] and, next establish strong monotonicity in Theorem HJ thereby concluding generic convergence for this 
system. 

Theorem 2>.IfQ is bipartite and Assumption\2\holds, then the system (0)-(0 is monotone. 

Proof: Let I c {1 , • • • , AO and I' = {1, • • • , AO - 1 be defined as in Proposition \T\ and suppose that in (O, the 
cells are indexed such that i - 1,- • • ,Nj belong to I, and i = Nj + 1, • • • ,N belong to I' as in the proof of 
Theorem^ where Nj is the cardinality of set I. Let X 1 := [jc iT • • • x N ^ ] T , X 1 ' := [x N ^ +lT ■ ■ ■ x nT ] t , and 
define U 1 , U 1 , Y 1 , Y 1 ' similarly. Then, the interconnection condition dU) becomes: 

U 1 = {Pn®hn)Y r (59) 
U r = (P2i®I m )Y J (60) 

where Pn and P 2 \ are as in (l37l) . A block diagram illustrating this interconnection is depicted in Figured] 



A' 1 -^A- 1 SiC 2 A- 1 
A' 1 





B 2 



- -A^C l A1 l B i C 2 A 2 l B 2 



(57) 



u 1 


x l 






















x 2 


Y 1 
=> 


P 2 i®I 


U 1 ' 






Y r 
=> 


P\2®I 




1 — > 








> 










X N Z 








x N 



















Figure 1 : A block diagram for the system (O-© when the contact graph is bipartite and the corresponding intercon- 
nection matrix P is decomposed as in d37b . 

To prove the monotonicity of this feedback system, we establish the monotonicity of the feedforward 
system with input H :- U 1 and output J/ := (.P\ 2 ®I m )Y r ': 
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Claim: The feedforward system in Figure \J} with input II and output J/ is monotone with respect to the 
positivity cones K u = K Y = R'^ J , and K x = K N i x {-K} N ~ N i. 

The theorem follows from this claim because a monotone input-output system, where the inputs and outputs 
are ordered with respect to the same positivity cone, is monotone when the output is connected to the input 
with unitary positive feedback H = J/ (see the first part of the proof of lfT2l Theorem 2]). 

To prove the claim above, we take two input signals satisfying ^(t) < f U(f) for all t > 0, which means 
that u\t) < u'(t), i e J, with respect to R™ Q . Likewise, we let X(0) < X(0) with respect to the cone K x - 
K Nl x {-K} N ~ Nl , which means that x'(0) < x'(0) for i e I and x'(0) > x'(0) for i e I' with respect to the cone 
K. It follows from Assumption |2] that: 

x\t)<x\t) Vt>0 iel. (61) 

Moreover, since x < x. implies h(x) < h(x) with respect to by Assumption |2l we conclude Y 1 (f) > Y T (t) 
with respect to R.> J . Because P21 is a nonnegative matrix, (l60l ) implies U r (t) > U r (t) which means that 
u'(t) > u'(t) for all t > 0, i e I'. As noted above, x'(0) > x'(0) for i e I' and, hence, another application of 
Assumption[2] yields: 

x\t)>x\t) Vf>0 iel'. (62) 

Since ([SB and ([62]) hold with respect to K, we conclude that X(t) < X(t) for all t > with respect to K x = 
K Nj x {-K} N ~ Nf . To conclude the proof of the claim, we need to show that X < X implies J/ < }f. Indeed, 
the former implies x' > x' for i e I' and, it follows from Assumption|2]that h(x') > h(x') with respect to R™ Q . 
Thus, Y 1 ' < Y 1 ' with respect to M.™^~ Nl ^ and, since Pu is a nonnegative matrix, we conclude J/ < }f with 
respect to R'"^ J . □ 

To establish strong monotonicity, we need additional excitability and transparency conditions, as defined 
inElEl: 

Definition 2. The monotone system x - f(x,u), y - h(x) is said to he excitable ifx(0) < x(0) and u(t) < u(t) 
for almost all t > imply x(t) <sc x{t) V? > 0. It is said to be transparent if u(t) < u(t) and x(0) < x(0) imply 
y(t) « y(t) Vt > 0. 

Since inputs and outputs are ordered with respect to orthants (K u = W" Q and K Y = -K u ) in Assumption 
121 here we give a less restrictive definition of excitability (transparency) which requires that this property 
hold with respect to a particular component of the input (output) vector: 

Definition 3. The monotone system x = f(x, u), y = h(x) is said to be excitable by the kth input ifx(0) < x(0), 
u(t) < u(t) and Uk(t) < Uk{t)for almost all t > imply x(t) <sc x(t) t > 0. It is said to be transparent from the 
kth output ifu(t) < u(t) and x(0) < x(0) imply yk(t) < yk{t) V? > 0. 

Assumption 3. There exists k e {1, • • • ,mj such that system (0) is excitable by the kth input and transparent 
from the kth output. 

Theorem 4. If in addition to the conditions of Theorem \3\ Assumption \3\ holds, then ©-(0 is strongly 
monotone. 

Proof: We need to show that X(0) < X(0) implies X(t) <sc X(t) for all t > with respect to the cone K x - 
K Nj x {-K} N ~ Nj for which monotonicity was proven in Theorem [3] By this monotonicity property, we 
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already know that X(t) < X(t) for all t > 0, and Assumption |2] implies y\t) > f{t) if i e I, y'{t) < y'(t) if iel', 
both with respect to Because P12 and P21 i n d59b-d60b are nonnegative matrices, we conclude: 

u\t) < u\t) i e I, u\t) > u\t) iel'. (63) 

Next, note that X(0) < X(0) means x'(0) t x'(0) for at least one i e {1, ■ • • ,N}, say i*. Thus, with respect to 
the cone K: 

/ (0) < f (0) if f e I, / (0) > f (0) if f el'. (64) 

Using (I63l-(l64l and the transparency assumption from the kth output, we conclude that the following holds 
for all f > with respect to the standard order induced by R>o: 

y'k (0 > f k (t) if i* 6 1, y'k (0 < fk « if ** £ r ■ (65) 

Now, pick an arbitrary r e {1, • ■ ■ ,N], and note from the connectedness of the contact graph Q that a path 
of adjacent cells />, I = 1, • • • ,M exists such that i\ = i* and im = ft. Since i% is a neighbor of /1 = i*, for all 
?>0, 

«J (0 > u l l (t) if i* e I, u h - (t) < fij (?) if f e J'. (66) 
Since is bipartite, i* e I means 12 £ X', and /* e J' means ^ £ Thus, from X(0) -< X(0): 

x fa (0) > x 1 ' 2 (0) if f e I, x h (0) < x h (0) if f e J' . (67) 

From the excitability assumption by the kth input, (l66l ) and (I67T ) imply: 

x ! ' 2 (0 » x i2 (/) if f g I, x h (t) « x ! ' 2 (0 if i* e I' (68) 

V? > and, from transparency, the following holds with respect to the standard order: 

< fl{t) if f e I, > if f e I'. (69) 

Continuing recursively, we conclude that (|66t-(f69T> hold for if, £ - 3, •• ■ ,M, with the inequalities reversed 
when £ is odd. In particular, (I68T ) becomes: 

(-l/x' v (f) » (-1)^(0 if f e J, (-l/x ,v (?) « (-l/x !> (0 if f e J'. (70) 

Since is bipartite, if M is even, /* -i\ el means ft - im £ and i* e I' means ft e I. Likewise, if M is 
odd, i* e I means ft e I, and i* e I' means ft e I' . Thus, d70l with ^ = M gives: 

/ (0 » x? (t) if j f e I' , / (?) «: / (?) if j f e J. (71) 

Since this inequality holds for each ft 6 {1,- • • ,N}, we conclude X(t) <sc X(t) as desired. □ 

In preparation for the examples in the next section, we now review a graphical test to ascertain excitabil- 
ity and transparency, given in lfT2l . Suppose the system x = f(x,u), y = h(x), xef c R", uef c R m , 
y e & c W, is such that, for each j + k, dfj(x,u)/dxf < is either identically zero, strictly positive, or strictly 
negative for all (x, u) e X x . Likewise, dfj(x,u)/duk and dh/ c (x,u)/dxj have the same sign definiteness 
property for each j and k. Associate to this system a directed incidence graph with vertices x\,--- ,x n , 
U\, • • • , u m , yi , • • • ,y p . A directed edge is drawn from x^ to xj, j ± k, if dfj(x, u) jdx^ is nonzero, from to xj 
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if dfj(x,u)/duk nonzero, and from xj to y^ if dhj c (x,u)/dxj is nonzero. The following lemma, adapted^! from 
lPT2l . proves excitability and transparency for systems that are monotone with respect an orthant cone: 

Lemma 6. Suppose the system x = f(x, u), y = h(x) is monotone with respect to an orthant cone and admits 
an incidence graph according to the rules described above. This system is excitable by the kth input if each 
state xj is reachable through a directed path from u^, and transparent from the kth output if a directed path 
exists from each state xj to yu- 

5 Examples 

5.1 A Class of Systems that Encompasses the Notch Signaling Model of (3]] 

As a special case of (f3]>, consider the single-input, single-output system: 

x\ = -y\x\+gi(x' 2 ) 

4 = -7/4 + ^(4-,) (72) 

X'n = -JnX l n + g„{u') 

y = A 

where, for / = 1, x'. > denotes the concentration of species j in cell i, y> > represents the corre- 
sponding degradation rate, and gj : R>o — > R>o is a continuously differentiable function. 

The reference |3] studied (1721 ) for n = 2 species, as a rough model for Notch signaling where the 
membrane-bound Delta ligands bind the Notch receptors in adjacent cells. This leads to the cleavage of 
Notch and the release of its intracellular domain which then serves as a co-transcription factor that inhibits 
the production of Delta in the same cell. Thus, in (1721 ). x\ represents the concentration of Delta and xj rep- 
resents the concentration of the co-transcription factor obtained from Notch. The function g\{-) is assumed 
to be decreasing since the co-transcription factor inhibits the production of Delta, and gi(f) lS assumed to be 
increasing since Delta activates the production of the co-transcription factor in adjacent cells. 

The reference O proved the emergence and stability of patterns for the case oiN = 2 cells, and observed 
the patterning behavior for N >2 by numerical simulations. A detailed bifurcation analysis is performed for 
this model in ifTTl . again for N = 2. In Proposition [2] below, we show that the results of the present paper 
are applicable to the model (1721 without restrictions on the number of species and cells. In particular, the 
instability criterion for the homogeneous steady-state in Theorem Q] makes use of the spectral properties 
of random walks and, unlike (3l HH which analyze this steady-state for specific arrays, is applicable to 
arbitrary graphs. Likewise, our study of checkerboard patterns in Proposition [Hand Theorem |2] generalizes 
the statements in (3l for two cells to bipartite graphs of arbitrary size. In addition, we establish monotonicity 
properties for bipartite graphs, thus revealing the global behavior of the solutions. 

2 Theorems 4 and 5 in 1121 give analogous tests for excitability and transparency with respect to Definitionf2] Theorem 4 requires 
that each state be reachable from each input through a directed path, and Theorem 5 stipulates that a directed path exist from each 
state to each output. The statement in Lemma|6]for transparency from the kth output follows directly from Theorem 5, by taking 
to be the only output of the system. The statement for excitability by the kth input follows from a straightforward modification of 
Theorem 4: Read the second part of the proof of Theorem 4 by replacing j* with k. 
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Proposition 2. System A72\) satisfies Assumption^ If an odd number of the functions gj(-), j = 1, • •• ,n, are 
nonincreasing and the rest are nondecreasing^ then it satisfies Assumption [2] as well. If in addition, the 
nondecreasing (nonincreasing) property is strengthened as: 

g'j(s)>0 (g'j(s)<0) Vs>0 j=\,--,n, (73) 

then Assumption\3\also holds. 

Proof: We first prove that Assumption Q] holds. Given u* > 0, the unique steady-state x* of (l72l) is given by: 



x n = 7n Sn{u ), Xj 

In particular, the input-output characteristic is: 



7j 8j(x* j+1 ), j = n- 1,— , 1. 



T(-):=y- x l gi(y- 2 l g2(---(y- x g n {-)))). 



(74) 



(75) 



The Jacobian matrix: 



df(x, u) 
dx 



-7i fe) 



-72 



Cl^") 
-7n 



(76) 



is upper-triangular with negative diagonal entries -yj, j ,n, and, hence, Hurwitz. This means that the 
determinant condition © holds and the steady-state x* is asymptotically stable. Note from (1721 that x n (t) 
exists for all f > and converges to y^gniu*). Applying a similar argument recursively for j — n - 1, •• • , 1, 
we conclude that x* is globally asymptotically stable. 

To show that Assumption |2] holds, we first select numbers e y - e {0, 1}, j = ,n, according to the 
following rule: Set e„ = if g n (-) is nondecreasing, and e n = 1 if g„(-) is nonincreasing. Then, for j = 
n - \,n -2, ■ •• , 1, set ej = ej+i if gj(-) is nondecreasing, and ej + ej+\ if gj(-) is nonincreasing. It follows 
from this construction that, Vs > 0, 



(-ir"g' n (s) > 0, (-ir^ +1 g'j(s) >0, j= ,n-l. 



(77) 



Since an odd number of the functions gj(-) are nonincreasing, the selection of the numbers ej above yields 
e\ - 1. Thus, an application of Lemma [5] with 8\ = and pi = I shows that the system (1721 ) is monotone 
with respect to K u = R> , K x = [x 6 R" | (-l)^xy > 0}, ^ y - R< , as in Assumption [2] 

To show that Assumption [3] holds, we apply the test in Lemma [6l The incidence graph for system d72l ) 
consists of the single path mhi„h x„_i h-> • • • h> xi i-> y, which means that any state is reachable from the 
input, and the output is reachable from any state. Thus, the system d72l is excitable and transparent. □ 



5.2 A Multi-Input, Multi-Output Model for Notch Signaling 

We now study the following system adapted^ from the lateral inhibition model in (H : 

N' = ^-yN i -kN i (D j ) i 
V = giS^-yD'-kD'iN^i 
S i = -yS l + kN i (D i ) i . 



(78) 
(79) 
(80) 



If one of the functions is constant, then one can count it as either nonincreasing or nondecreasing. However, this situation is of 
no interest in this paper, since the input-output characteristic l !75t is constant and, therefore, Theorems Q] and [2] are not applicable. 
4 The equation corresponding to J801 > in (5) includes a Hill function of N'(D-i}j instead of the linear term used here. 
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Here, A? 7 > 0, D' > 0, S ' > are the concentrations in cell i of the Notch receptor, Delta ligand, and a signaling 
protein activated by the binding on Delta and Notch, k > 0, y > 0, /? > 0, g : R>o — > R>o is continuously 
differentiable, and decreasing since the production of Delta is inhibited by the signaling protein. The notation 
O; denotes the average of the quantity within brackets over all cells adjacent to i. Unlike the model of Q 
discussed in the previous subsection, (l78l) incorporates the Notch receptor. 
We let: 

u\ := {D j )i, 4 := (N j )i, x\ = N\ x\ = D', x\ =N' + S \ (81) 

and rewrite {78]>-([80]> as: 



l 3 
.V', 



fi-yx\ -kx\u\ 

g(x' 3 — jcj) — yx l 2 - kx l 2 u' 2 

-J A +fi 



which is of the form © with X = {x € R 3 | xi > 0,x 2 > 0,x 3 > xi], 9/ = <¥ = 



I 2 

So- 



(82) 
(83) 
(84) 
(85) 
(86) 



Proposition 3. The system fl82D-fl<S6D, where k > 0, y > 0, yS > 0, and g : R>o — > R>o is continuously differ- 
entiable, satisfies Assumption^ Ifg(-) is nonincreasing, then it also satisfies Assumption^ Ifg'(s) > Ofor 
all s>0, then Assumption\3\holds for solutions in the forward invariant subset of & where x\ > 0, x 2 > 0, 



Proof: Given u* >0,u* 2 > 0, the unique steady-state of ([82])-([84]) is given by: 



x, , 

y + ku*. 



4-"-, 

7 



g(*3~*P 

y + ku* 



(87) 



and the Jacobian matrix: 



df(x, u) 



dx 



(x,u)=(x* ,u*) 



-J. 

-g (x 3 ■ 





-ku*. 



x*) 





-y - ku*, 






?'(^-xp 

-y 



(88) 



has the negative eigenvalues -y-ku* v -y-ku* 2 , -y, and is thus Hurwitz. It is clear from (1821 and (1841) 
that x\(t) and x 3 (i) converge to x\ and x*. This means that the first term in ([831 converges to g(x* 3 -x*), 
from which we conclude that x%(i) converges to x*. Thus, x* is globally asymptotically stable and all other 
statements in Assumption Q] hold. 

To verify Assumption |2l we note that: 



df\ , ^ n df 2 df 2 , df 2 

— = -kxi < 0, — = -kx 2 < 0, — = -g (x 3 - xi) > 0, — 

OU\ OU 2 OX\ OX3 



dh 1 dh 2 

? '(X3-X!)<0, -± = 1,-1 = 1. 

OX 2 OX\ 



(89) 



Thus, Lemma[5]holds with 5\ = S 2 = 0, pi = p 2 = 1, e\ = e 2 = 1, e 3 = and, thus, we conclude monotonicity 
with respect to the orthants K u = R^ , K Y = -K u , K x - {x e R 3 \xi < 0,x 2 < 0,x 3 > 0}. 

To show that Assumption [3] holds, we apply the test in Lemma [6] The incidence graph for the system 
(I82l-(l86l restricted to the subset of X where x\ > 0, x 2 > 0, X3 = x* 3 is as in Figure [2] From Lemma [6l 
we conclude that the system is excitable by u\, since a directed path connects u\ to both x\ and x%, and 
transparent from y\ , since a directed path connects both x\ and x 2 to y 1 . □ 
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Figure 2: The incidence graph for system (l82li-(l86ll. constructed as in Lemma|6] 

Note that the restriction x 3 - x* allowed us to drop xj, which is not excitable by either input, from the 
incidence graph in Figure |2] Likewise, the restriction Xj > 0, j - 1,2, is critical for the sign-defmiteness of 
dfj/duj = -kxj, which made it possible to direct an edge from Uj to Xj. Because the subset of 3£ denned 
by these restrictions is forward invariant and attractive, the w-limit sets of all solutions starting in 3£ lie in 
this subset. Thus, strong monotonicity on this subset, established by Theorem H] when the contact graph is 
bipartite, allows us to conclude generic convergence on SC . 

We emphasize that the assumption of identical degradation rates for N and S in (I78l)-(l80l) is essential 
for the change of coordinates that lead to (I82l-(l86l) and that allowed us to conclude monotonicity using 
Lemma[5]for orthant orders. It would be interesting to investigate whether monotonicity can be established 
for nonidentical degradation rates. 
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